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PRLIACfc 


As  military  aircralt  lly  with  more  and  more  stores,  the  problem  ol  predicting  the 
critical  flutter  speed  becomes  more  and  more  complicated.  Apart  Irom  the  lack  ol 
aerodynamic  knowledge  which  makes  it  almost  impossible  to  predict  the  unsteady 
aerody  namic  forces,  the  non-linear  behaviour  of  wing-stores  configurations  leads  to  fresh 
difficulties  m the  interpretation  ol  ground  vibration  tests  and  (light  flutter  tests. 

Mr  llaidl's  paper  helps  in  the  understanding  of  these  difficulties  with  the  support  ol 
ground  and  llight  experience  on  modern  aircraft.  Ilis  contribution  will  be  ol  great  value 
lor  all  concerned  NATO  countries. 


<>.(  OUPRY 

Chairman.  Sub-Committee  on 
Aeroelasticity  and  Unsteady 
Aerodynamics 
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SUMMARY 


I- samples  of  non-linear  vibration  behaviour  in  ground  resonance  tests  of  an  aircralt 
are  shown  Model  tests  lor  a simplified  system  with  non-linear  properties  have  been 
performed  to  study  the  effects  ot  friction  and  backlash  with  respect  to  ground  resonance 
test  and  flight  flutter  test 

With  symmetric  and  asymmetric  non-linear  stillness  characteristics  effects  ol 
amplitude  dependent  frequencies,  mode  coupling,  mode  asymmetries  and  the  consequences 
in  parameter  identification  in  vibration  tests  are  pointed  out  and  discussed. 

In  case  of  flutter  critical  modes  the  problems  of  apparent  damping  caused  by  non- 
linear system  properties  are  shown  and  recommendations  are  given  to  reach  a represen- 
tative flutter  clearance  with  respect  to  this  non  linear  system  behaviour. 
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PHOHU.M  SbHVEY 

lhere  have  been  many  improvements  in  equipment  and  test  methods  for  ground  resonance 
tests  and  flight  flutter  tests.  The  determination  of  the  dynamic  characteristic  of  real 
aircraft  structures  is  still  difficult  in  case  of  non  linear  system  behaviour. 

Some  examples  taken  from  previous  tests  and  investigations  may  help  to  give  a problem 
survey. 

The  first  example  refers  to  a ground  resonance  test  on  a F 104  G MS tarf i ght er"  air- 
craft, described  in  AGAR1)  Keport  No.  592,  see  fhef.  l]  . Two  mul  t mat  ionale  comparative 
programs  have  been  performed  to  define  and  improve  the  aerodynamic  and  e las tomechanical 
basis  for  the  evaluation  of  aircraft  aeroelastic  stability.  Test  results  of  intensive 
ground  resonance  tests  indicate  the  common  trend  of  slightly  reducing  frequency  with 
increasing  amplitude  (fKef.  l],  Table  23  and  2**).  A strong  non  linear  behaviour  has 
been  found  for  a fin  torsion,  tailplane  coplanar  rotation  mode  (r  = h ) . With  increasing 
amplitude  the  resonance  frequency  drops  from  12.7  Hz  to  a minimum  of  10.2  Hz  but  in- 
creases to  11.2  Hz  with  higher  amplitudes.  Such  non  linear  effects  may  be  caused  by 
backlash,  friction,  complex  actuator  stiffness,  liquid  loads,  or  engine  installation. 

An  investigation  of  dry  friction  effects  to  an  aileron  rotation  mode  of  a glider 
has  been  performed  by  O.N.E.R.A.  [ He f . 2]. 
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Fig.  1 Influence  of  an  Auxiliary  Excitation  (60  Hz)  on  the  Aileron  Rotation  Mode 
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Fig.  4 Effects  of  Non  Linear  Stiffness  to  the  Wing  Aileron  Flutter  Stability 

f ig,  '*  reVwals  the  mast  critical  behaviour  of  the  system  with  the  shown  degressive 
stiffness  characteristic.  This  flutter  behaviour  is  very  dangerous  because  there  is 
no  amplitude  limitation  after  the  system  reaches  the  initial  amplitudes  of  the  shown 
stability  boundary. 

All  these  examples  may  illustrate  the  problems  in  the  ground  resonance  test  to  de- 
fine the  amplitude  dependent  dynamic  characteristics  of  a structure  and  if  possible 
to  lind  a representative  "linearisation”  at  a reference  point  close  to  operational 
cond 1 t i ons • 

I he  following  examples  from  a recent  ground  resonance  test  on  an  aircraft  with 
sweepable  wing  and  underwing  stores  may  show  the  limits  of  this  linearisation  technique 
and  possible  consequences  to  the  test  philosophy. 

I want  to  acknowledge  the  work  of  Mr.  John  B.Cox  from  Mechanical  Test  Department 
of  MAC-Warton,  who  has  contributed  considerably  to  the  A/C  test  examples  which  are 
shown  in  this  paper. 


Non  linear  Effects  in  A/C  Ground  Resonance  Test 

The  non  linearities  which  will  Vie  dicussed  here  are  effects  of  backlash  and  friction, 
influencing  the  dynamic  system  behaviour.  Fig.  5 gives  a sketch  of  the  simplified  vi- 
bration system. 
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Fig.  5 Sketch  of  the  Simplified  Vibration  System 
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I he  wings  are  sweepable  in  teflon  coated  bearings  and  are  driven  by  wing  sweep 
actuators,  1 he  friction  moment  in  this  teflon  bearings  is  relatively  high  (^2  kNm) 
because  of  "shrunk  in  procedures"  and  static  loads,  which  is  much  higher  than  usual 
available  excitation  moments  with  tiifl  equipment.  As  long  as  the  external  exciting 
moments  or  t fie  internal  moments  of  the  vibrating  system  in  wing-yaw  direction  are 
below  the  friction  moment,  actuator  backlash  and  actuator  spring  stiffness  have  no 
u»  limited  influence.  1 he  system  properties  are  changed  when  the  bearing  is  moving. 

I he  wing  slot  sealing  is  introducing  additional  stillness  and  friction,  varying 
with  wing  sweep  angle. 

The  attachment  ol  the  actuator  at  the  wing  carry  through  structure  has  a lateral 
deflection  characteristic  which  lias  an  influence  on  asymmetrical  modes  whereas  1 or 
symmetrical  modes  the  lateral  forces  normally  are  compensated. 

Considering  t fie  store  yaw  mode,  there  is  a similar  friction  coupling  mechanism  in 
the  teflon  bearing  of  tfie  pylon  spigot.  This  friction  and  backlash  in  the  pylon  sweep 
drive  rod  mechanism  generates  similar  non  linear  system  behaviour  as  discussed  before. 
I fie  attachment  of  the  store  to  tfie  pylon  creates  additional  backlash  and  friction 
e t t e c t s . 


Fig.  6 Effect  of  Input  Power  on  Resonance  Frequency  for  Wing  Yaw-Mode 


Fig.  6 shows  measured  decrease  of  the  resonance  frequencies  with  forcing  level  for  the 
wing  yaw  mode  (curve  A).  After  a back  and  forward  sweep  of  the  wing  a change  in  fre- 
quencies have  been  found  (curve  B ) . 


This  trend  of  frequency  decrease  is  given  by  increasing  effects  of  short  movements 
in  the  bearing,  changing  the  stiffness  and  damping  of  the  system.  It  should  be  mentioned 
fiere , that  the  friction  coefficient  for  teflon  in  static  and  gliding  condition  is  about 
the  same.  That  means  that  the  part  time  movements  during  one  vibration  cycle  are  limited 
for  durating  system  moments  higher  than  the  friction  moment.  By  this  reason  an  auxiliary 
excitation  as  mentioned  before  cannot  "break"  this  friction. 

Effects  of  a preload  in  wing  aft  direction  to  the  wing  yaw  mode  can  be  seen  in  the 
time  histories  of  the  actuator  forces,  Fig.  7a. 


Fig.  7 Time  History  of  the  Actuator  Force 

Wing  Yaw  Mode  Excited  with  90  N on  Outboard  Pylon  Station 

With  compensated  preload  (wing  suspended)  the  resonance  amplitudes  with  nearl\  the 
-*ame  excitation  moment  are  reduced  as  shown  in  Fig.  7 h.  This  non  expected  behavioui  c 
be  explained  by  reduced  gliding  effects  with  preload. 

Ihe  non  symmetric  stiffnes  characteristic  is  indicated  in  the  wave  form  of  the 
port  actuator  force  in  l lg.  7a.  The  energy  consuming  relativ  movements  in  the  wing 
pivot  bearing  (friction  coupling)  therefore  seem  to  be  limited  to  the  wing  forward 
oscillation  cycle. 

In  the  symmetric  case  without  preload  there  are  gliding  effects  in  both  directions 
causing  higher  system  damping  and  reduced  amplitudes,  see  Fig.  7h . 

More  severe  effects  of  friction  and  backlash  have  been  found  for  the  store  modes. 
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Fig.  8 Vibration  in  Frequency  of  Low  Frequency  Modes 


em  . 

'ara  lie 
la t 1 ng 
1 on  to 
of  t h 


he  mod 
fitted 
be  twee; 
eren  t 


iine  h 
ch  i n.vt 


7 


Is 

ww 

ExCtlAHUN 

\ * \ * \ !\  / ■ 

VELOCITY  Vy 

forward  jtore 

A A A A 


veloc.ty v, 

REAR 


f*\  .***  /*A 

J \*i  vy  V/ 

C SPLACfcVE  NT 
DAMPS*  PA Hi  1 


..  MEM 
DAMP)  » pah:  2 


GUD-'*  - TH  GN 


A . v 

DAMPER  PAHT  1 AND  2 


, .30-10  ' ' 

■ 0 ,Nrn) 

L-30O'2 

pZA 


'//.V.V.V.V.V.^v/.v.v/.v  )•  o iu 

-2A 


.'.\7.7A\7.\ 


. •7-1'  •■htf'r.,'/'  ,Vl£.i 


0 l6m/l 

-2A 

• UK)'3 

0 |»«dl 
-12  »'3 
.42 10'4 
0 i'«*» 

-w‘ 

. 1J  IS'3 
0 ;■«. 
-v,*)-3 


.«I0'2 

0 iNmi 

-6JD-10-2 

r6 0 

0 |cm/»J 

-60 

.60 
0 I cm/.) 
-60 

• 29B"3 
0 i '"I 
-29KT3 

i-4 


r* 


4210s 


-42B"' 


•2610"J 
0 i«ii 
-26  B'3 


Fig.  10  Time  Histories 


Friction  Moment  7,6  10  * Nm 
Excitation.  1,37  10  ^ Nm 

Resonance  Frequency:  3,80  Hz 


Fig.  11  Time  Histories 


Friction  Moment:  7,6  10'^  Nm 

Excitation:  3,67  10’*  Nm 

Resonance  Frequency:  3,62  Hz 


Ihe  trace  oi  the  relative  motion  in  Fig.  10  demonstrates  the  part  time  movements  in 
the  friction  coupling  during  a vibration  cycle.  As  long  as  there  is  a relative  motion, the 
decay  curves  show  the  effects  of  gliding  friction  together  with  structural  damping. 

After'  the  relative  motion  is  stopped  the  decay'  curve  shows  the  structural  damping 
only  • 

Fig.  11  and  Fig.  12  illustrate  the  change  in  the  decay  curves  for  increasing  exci- 
tation moments  and  the  same  friction  moment. 
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Friction  Moment:  7,6  10'zNm 

Fig.  12  Time  Histories  Excitation  6,1110 2 Nm 

Resonance  Frequency.  3,54  Hr 


Fig.  13  Time  Histories 


Friction  Moment:  1,6  10*' Nm 

Excitation:  7.34  10  ^Nm 

Resonance  Frequency:  3,72  Hr 
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1 he  t line  hi  -story  in  I 111.  13  show  s n c a to*  with  higher  1'riction  moment,  Similar  meas- 

urements with  backlash  in  the  spring  attachment  of  damper  part  1 indicated  a relatively 
rapid  stop  of  the  relative  movement  in  a position  somewhere  within  the  backlash  range. 

In  1 lg,  1 the  amplitude  dependent  effects  of  a constant  friction  moment  to  the 
apparent  system  damping  are  shown. 


Fig.  14  Apparent  Damping  vs  Exciter  Moment 
(g  = 1 ,3  % Structural  Damping  Included) 

With  increasing  amplitudes  the  clamping  reaches  a maximum  and  decreases.  This  behav- 
iour demonstrates  the  misleading  effects  upon  the  system  damping  expressed  in  % g and 
it  demonstrates  the  reason,  why  friction  damper  are  not  allowed  to  suppress  flutter. 

An  analytical  approach  of  this  system  behaviour  with  the  so  called  method  of  "Har- 
monic Linearisation"  was  performed.  An  example  of  the  calculated  resonance  amplitudes 
and  frequencies  compared  to  test  results  are  shown  In  Fig.  15. 


Fig.  1 5 Resonance  Amplitude  and  -Frequency  vs  Excitermoment 


The  actual  amount  of  backlash  and  friction  of  auch  a system  can  be  determined  by 
quasi  static  measurement  of  the  "hysteresis  diagram",  see  Fig.  l6. 


Fig.  16  Hysteresis  Diagram 


This  stiffness  characteristic  is  measured  with  an  exciter  moment  acting  on  the  store 
(frequency  0.02  Hz).  The  diagram  of  store  rotation  (curve  A)  and  rotation  of  damper  part 
1 (curve  B)  reveals  the  amount  of  backlash  and  friction  of  the  actual  system. 

Kltects  of  Friction  and  Backlash  to  Data  Analysis  of  Flight  Flutter  Tests 

Different  excitations  have  been  used  in  the  model  test  in  order  to  study  the  effects 
to  the  data  analysis  with  different  methods  [Ref.  5i  6j  . 
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Fig.  17  Time  History  of  a Test  with  Random  Excitation 
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!•  i g . 17  is  a time  history  ol  a test  with  random  excitation  which  shows  tin*  random 

excitation  signal,  the  lateral  response  Vy  at  forward  and  real  stole  and  the  relative 
movement  m the  friction  coupling.  The  model  has  a backlash  of  0,  1 mm  ( I , fi  • lO”-5  rad/ 

in  the  attachment  of  the  "pylon  sweep  rod"  springs. 
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Fig.  18  Time  History  of  a Test  with  Random  Excitation 


Fig.  18  shows  the  increase  of  relative  movements  in  the  friction  coupling  with  in- 
creasing excitation,  affecting  the  system  damping  and  resonance  frequency. 
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Fig.  19  Time  History  of  a Test  with  Sweep  Excitation 

Fig.  19  illustrates  the  dynamic  response  of  the  system  to  a frequency  sweep  input. 
Relative  motions  in  the  friction  coupling  appear,  when  the  excitation  frequency  is 
running  through  the  system  resonance. 
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Fig.  20  Time  History  of  a Test  with  Sweep  Excitation 

Fig.  20  shows  the  dynamic  response  of  the  same  system  excited  with  higher  force 
level.  There  are  considerable  time  variant  effects  to  the  system  damping  and  resonance 
frequency • 

In  a flight  flutter  test  with  artificial  frequency  sweep  excitation  and  superimposed 
random  excitation  there  is  a combination  of  the  shown  behaviour  for  random  and  sweep 
excitat ion. 


A mix  of  different  effects  of  the  friction  moment  and  different  resonance  frequen- 
cies in  the  response  signal  are  the  consequence  . 

Digital  data  analysis  with  statistical  methods  is  based  on  the  assumption  of  a linear 
system.  Changing  resonance  frequencies  with  amplitude  will  produce  an  apparent  damping, 
which  is  not  existent  in  the  system.  As  shown  before  in  Fig.  1*4, a friction  moment  is  in- 
troducing similar  misleading  effects  of  apparent  damping. 

Fig.  21  is  an  example  of  data  analysis  of  a random  excited  model  test  with  "backlash". 
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Fig.  21  Frequency  and  Damping  Analysis  of  Random  Excited  Model  Data 
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It  shows  the  computed  complex  plot  of  the  transfer  function  II  ( i <J  ) and  after  a 
"las*  Fourier  transform"  the  impuls  function  li^j  "nil  the  "log.  magnitude"  of  the 
unpulsfunct ton. 
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In  the  unpulsf unc t ion  there  is  no  more  an  indication  of  amplitude  dependent  fric- 
tion effects  and  frequency  change,  ihe  analyzed  apparent  danij>mg  is  considerably 
higher  than  the  .structural  system  damping. 
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ihe  shown  effects  of  friction  and  backlash  to  an  A/C  wing  yaw  mode  and  store  yaw 
mode  may  be  rather  extreme  examples  of  non  linear  effects,  hut  the  intent  oi  the 
studies  is  to  provide  information  leading  to  a better  understanding  of  non  linear 
system  behaviour  and  to  show  consequences  in  the  identification  ol  important  system 
p arame t ers . 

Considering  the  stiffness  characteristic  of  the  non  linear  "store  yaw  model"  as 
shown  before  in  I ig.  10  the  boundary  conditions  of  the  system  are  defined  by  pylon 
"fixed  by  friction"  for  small  amplitudes  and  pylon  spring  in  series  with  the  spring 
"sweep  drive  rod"  for  large  amplitudes. 


Between  these  small  and  large  amplitude  boundary  conditions  a decreasing  and  in- 
creasing effective  stiffness  with  amplitude  is  caused  by  the  amount  of  friction  and 
backlash.  The  minimum  of  the  effective  stiffness  is  a further  important  condition  in 
respect  to  parameter  variation  in  the  "linear"  calculation.  This  minimum  effective 
stiffness  is  strongly  dependent  on  change  in  friction  and  backlash.  Without  friction 
moment,  the  effective  stiffness  within  the  backlash  range  would  be  zero.  In  order  to 
have  a representative  variation,  differences  within  the  production  line  and  wear 
effects  must  be  covered. 


Asymmetries  in  the  stiffness  characteristic  caused  by  preload  give  a change  in 
effective  stiffness  compared  to  the  symmetric  case. 


Vibration  measurements  somewhere  within  the  backlash  range  cannot  be  used  to  assess 
the  range  of  necessary  parameter  variation  in  the  vibration  calculation  or  to  verify 
the  calculation.  In  such  a case  the  separate  determination  of  the  hysteresis  diagram 
is  a valuable  step. 


The  Aircraft  Ground  Resonance  Test  in  the  store  yaw  mode  and  wing  yaw  mode  has 
shown,  that  a test  approach  to  the  boundary  condition  "large  amplitude"  as  discussed 
before,  could  not  be  reached. 

The  model  tests  with  harmonic  excitation,  with  random  excitation  and  with  frequency 
sweep  excitation  can  provide  necessary  informatio.is  on  the  behaviour  of  the  non  linear 
system • 

The  effects  of  these  non  linear  system  properties  to  test  methods  used  in  ground 
resonance  tests  and  the  data  evaluation  in  flight  flutter  tests  have  been  demonstrated 
in  the  model  tests. 


As  shown  in  the  different  examples,  the  non  existent  apparent  damping  in  the  data 
analysis  caused  by  amplitude  dependent  resonance  frequencies  and  misleading  effects 
friction  depicts  a very  serious  problem  in  the  flight  flutter  test. 

Asymmetries  in  the  modes,  preload  effects  and  crossing  frequencies  of  the  store 
yaw  mode  with  store  pitch  mode  are  causing  additional  problems  in  mode  separation  and 
damping • 
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An  essential  step  with  respect  to  Mutter  stability  or  fatigue  is  to  determine 
whether  or  not  these  non  linear  features  are  important  in  the  flutter  behaviour  of 
the  A/C.  Following  aspects  have  to  be  considered: 

• modes  relevant  for  flutter 

. modes  coupling  with  relevant  flutter  modes 

• asymmetries  in  modes  affecting  the  flutter  speed 

• effects  of  apparent  damping  in  the  data  analysis 

Hetonunendations  to  achieve  a representative  flutter  clearance  in  case  of  severe  non 
linearities  are: 

• careful  parameter  identification  in  static  tests 

• calculations  with  corresponding  parameter  variations 

• wind  tunnel  tests  to  find  out  the  influence  of  asymmetries  and  non  linear 
features  to  flutter 

• non  linear  analytical  investigations 

• ground  resonance  tests  and  flight  flutter  tests  with  a modified  test  airplane 
(reduced  damping  and  backlash)  to  cover  the  worst  condition. 
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